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INTRODUCTION 


The self-induced on-orbit molecular contamination around the Space Station 
(SS) is of interest to anyone (user) anticipating doing experiments on or in 
the vicinity of the SS. Aerospace engineers need to design a SS propulsion 
system that keeps the SS in a stable orbit and at the same time does not allow 
the propellant gases to interfere' with the experiments of the user. Another 
engineering consideration that has to be addressed in the SS design is a 
mechanism for venting waste gases. 

One scenario that might accomplish the above requirements is to use an 
electrothermal propulsion system, resistojet, that will thrust continuously in 
the hundreds of mi 1 1 i-newton range which will provide a constant altitude for 
the SS with a low g environment. As a first attempt to understand the 
contamination from such a propulsion system, a point source model was 
developed in the summer of 1985 by Riley ^ at NASA - Lewis Research 
Center. This paper gives the numerical results of the point source model. 
Number column densities for C0 2 are presented in Table 1 as a function of 
direction of observation (line of sight), temperature of the exit gas, and 
mean exit velocity. One set of number column densities for N 2 is shown in 
Table 1 to contrast the effect of molecular weight. All the results are for a 
constant exhaust rate of 5,000 kg/year. 


THE MODEL 


Propellant gases from the jet nozzle are assumed to come from the jet exit 

plane with a mean macroscopic velocity*? relative to the nozzle. The gases 
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are assumed to be exhausted into a vacuum. Thus, there are no collisions with 
ambient molecules. At a distance of a few nozzle exit diameters from the 
nozzle, the gases will appear to come from a point. 

The mathematical model developed ( using kinetic theory of gases ) allows 
one to calculate at any position and time the number density around a point 
source that has a Maxwellian velocity distribution. 


Using the Boltzmann Equation one can write (see appendix A for 
nomenclature). 

it + TT - v r f = S(r) Q(t)£ <■*) 

( 1 ) 

where 






2fF4(T , 



( 2 ) 


and 


*0 u) ^ 

( 3 ) 

The assumptions in writing down equation (1) above are: (a) no external 


forces are exerted on the molecules, (b) the self-interaction terms are 
neglected, and (c) the source expands into a vacuum. . 


Using a constant flux source Q<t)=Q 0 and the Fourier-Laplace integral 


transform method, the result for the number density is 
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where is defined by 

r*/U 

where r is the location of the point that the number density n is to be 
calculated and "u is the mean exit velocity. Both quantities are relative to 
the point source nozzle. 


CALCULATIONS 


In order to make the number density calculations as a function of 0, 


defined above, the erfc function in equation (4) was approximated (2) by 

Mf feto * ~£i -+-*371 

v 13 y'S' 


73to 'jfico 

for small values of the erfc argument, and by 
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for large values of the erfc argument. The transfer from equation 5 to 
equation 6 was made when x was greater than 1.68. 


The point of observation was assumed to be 50 meters directly above the 
point source (see figure 1). Once one knows n(r,t), the number density, it 
is a simple matter to calculate the number column density (NCED by integrating 
numerically along a particular line of sight (LOS). 

The number column densities as a function of the LOS were calculated for 
different temperatures at the exit plan (temperature of the gas at the point 
source) and different magnitudes of the mean exit velocity "if. The 
calculations were performed for the steady state situation for continuous 
firing. Therefore, the r/ut terms in equation (4) were assumed to be zero. 

The number of C0 2 molecules per second, Q q , exiting from the nozzle 

21 21 

was assumed to be 2.2x10 . For N 0 Q =3.5x10 . These values were 
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obtained by assuming Q q = 5,000 kg/year. 

RESULTS 


Table 1 gives the results for the NCD calculations as would be seen by an 
observer located 50 meters above the point source as indicated in fig. 1. The 
observation angle, f>, is defined in fig. 1 to be the angle between the 'line 
joining the observation point and the source and the line of sight. (For 
example, an angle of 270° for p would mean a line of sight in the direction 
parallel to "u . ) 

The results of the NCD calculations are given in table 1 as a function of 
the temperature of the source, the mean exit velocity, and the line of sight. 
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All calculations were done in the plane containing a line connecting the 
source and observation point and the mean exit velocity vector. 

Figures 2, 3, 5, 5, 6, and 7 respectively are graphical representations of 
the NCD's as a function of observation angle for the data presented as Case 1, 
2, 3, 4, 5, and 6 respectively in Table 1. The exit temperatures are the same 
for figures 2 and 3 with the mean exit velocity being 4,000m/sec. and 
l,000m/sec. respectively. Upon comparing the NCD's as a function of LOS in 
figures 2 and 3, one observes that higher column densities are prevalent at 
larger values of the angle between the LOS and the mean exit velocity vector 
which is at 270 0 in figures 2 and 3. This implies that high mean exit 
velocities are necessary to increase the acceptable viewing solid angle. The 
mean exit velocities are the same for figures 3, 4, and 5 with the exit 
temperatures being 1600°K, 800°K, and 100°K respectively. 

When one compares figures 3, 4, and 5, the effect of the exit temperatures 

on the NCD's is clearly seen. Low exit temperatures tend to keep the exit 

plume from expanding in a direction perpendicular to the mean exit velocity. 

This is clear from the fact that the NCD's are larger for observation angles 

above 270° (which are for lines of sight through the exit plume) for the 

lower temperature cases. This implies that low exit temperatures are 

important in order to increase the acceptable viewing solid angle around the 

1 5 2 

SS. If one assumes the acceptable value of 1x10 / m for NCD's, then 

figures 6 and 7 show that any LOS would have unacceptable values for NCD's for 
a zero mean exit velocity. These zero mean velocity cases could represent a 
point source leak. 
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To gain better insight into how the number density function at different 
positions around the nozzle contributes to the LOS number column densities, 
equidensity contours as a function of the angle are given in f-igures 8, 9, and 
10 for an exit temperature of 1600°K and a mean exit velocity of 4,000m/sec. 
(Case 1 in Table 1). Figures 11, 12, and 13 are similar plots except the mean 
exit velocity is l,000m/sec. (Case 2 in Table 1). Figures 14, 15, and 16 are 
also a similar set except for a mean exit velocity of l,000m/sec. and a 
temperature of 100°K (Case 4 in Table 1). Figures 17, 18, and 19 are 
equidensity plots for zero mean exit velocity (Case 5 in Table 1). 

Figures 20, 21, and 22 are presented to show how the equidensity contours 
are elongated (for a fixed mean exit velocity) as the temperature approaches 
zero. All three plots show how the gas expands in a direction perpendicular 
to the mean exit velocity (positive x-axis). Notice the maximum distance (y 
value) occurs at larger values of x for the lower exit temperatures. These 
plots are not to imply that C0 2 will remain a gas at these low temperatures 
as it exits from a real nozzle, but suggests how an ideal gas plume will 
behave as the temperature approaches zero with this model . 

DISCUSSION 

If the allowable limit for the NCD's for C0 2 is assumed to be IxlO 15 
molecules/m, then any observation angle greater than 270° (parallel to 
the mean exit velocity vector - ?) would be unacceptable for all the cases shown 
in Table 1. A study of Table 1 and the figures presented reveal, as one might 
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expect from the kinetic theory point of view, the need to have high mean exit 
velocities at low exit temperatures to confine the plume in such a manner as 
to maximize the allowable viewing solid angle. 

For a mean exit velocity of 2000m/sec. and a temperature of 100°K, a 
maximum observation angle limit of 260° is suggested by the data displayed 
in figure 23. For lines of sight at an observation angle greater than 260°, 
the NCD's will be greater than lxlO^/m^. 


TABLE 1 




ANGLE COLUMN DENSITIES 

(deq) ( part i c 1 es/m A 2 ) 

4 >~ 



CASE 1 

CASE 2 

CASE 3 

0 

30 

5. 7E+03 

1.6E+15 

3.4E+14 

60 

3. 0E+03 

7. 3E+14 

1.6E+14 

90 

2.4E+03 

5.4E+14 

1.2E+14 

120 

2.6E+03 

5. 1E+14 

1. 1E+14 

150 

4.0E+03 

6. 1E+14 

1.4E+14 

180 

1 . 2E+04 

8.5E+14 

2. 3E+14 

210 

6. 6E+06 

2.0E+15 

7.8E+14 

240 

3.8.E+12 

5. 7E+15 

4.7E+15 

270 

1.5E+16 

1.5E+16 

2. 1E+16 

300 

3.6E+16 

2.9E+16 

4.6E+16 

330 

360 

6.3E+16 

5.9E+16 

8.8E+16 



CASE 4 

CASE 5 

CASE 6 

0 

30 

2.3E+04 

6. 1E+16 

8.0E+16 

60 

1.2E+04 

2.SE+16 

3.5E+16 

90 

9.8E+03 

1.8E+16 

2.3E+16 

120 

1.1E+04 

1.4E+16 

1. 7E+16 

150 

1 . 6E+04 

1.2E+16 

1.5E+16 

180 

4.9E+04 

1.0E+16 

1.3E+16 

210 

2.7E+07 

1.2E+16 

1.5E+16 

240 

1.5E+13 

1. 4E+16 

1.7E+16 

270 

5.8E+16 

1 . 8E+ 1 6 

2.3E+16 

300 

1.5E+17 

2.SE+16 

3. 5E+16 

330 

2.6E+17 

6. 1E+16 

8.0E+16 

360 







CASE 

1: 

( CO 2 ) , TEMPERATURE 

CASE 

O * 

( C02 ) , TEMPERATURE 

CASE 

•C* ■ 

( C02 ) , TEMPERATURE 

CASE 

4: 

( C02 ) , TEMPERATURE 

CASE 

5: 

( C02 ) , TEMPERATURE 

CASE 

6: 

( N2 ) , TEMPERATURE 


1600 K, VELOCITY 4000 M/S. 
1600 K, VELOCITY 1000 M/S. 
800 K, VELOCITY 1000 M/S. 
100 K, VELOCITY 1000 M/S.. 
300 K, VELOCITY 0 M/S. 
300 K, VELOCITY T O -li/S. 
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FIGURE 7 

N2: TEMF-300 K, VELOCITY=Q M/5 
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FIGURE 8 

EQUIDENSITY PLOT ( 1 0"' 1 4/m "3) 


C02; TEWP= 1 6w K, VELOCITY =4C-C>0 M/5 
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FIGURE 9 


EQ U I DENSITY PLOT ( 1 0'- 1 5/m ~ 3 ) 


C02; TEMP= 1 SOO K, VELQC ITT =4QCO M/5 
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FIGURE 10 
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FIGURE 11 


EQ U I DENS ITT P LOT ( 1 1 4/m ~ 3 ) 


C02; 7EMP= 1 SCO K, VEL0CUT'=1 CvO M/5 
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FIGURE 12 


EQUIDENSITY PLOT (1tri5/m''3) 

C02; TEMP=15W K, VELOCITY = 1 WO M/5 
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FIGURE 13 

EQUI DENSITY PLOT (10~16/rrr3) 


TEMP=16W K, VELCCirr=1QOQ M/3 






FIGURE 14 

EQU I DENSITY PLOT (10^ 14/m 


C02; TEMP= 1<X> K, VELOCfTY= 1 000 M/' 
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FIGURE 15 


EQUI DENSITY PLOT (10''15/m^3) 


C02; TEMP=100 K, VELOC nT= 1 0OO M/5 
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FIGURE 16 


EQUI DENSITY PLOT ( ICT 1 6/m''3) 


C02; TEMP=KX> K, VELOCrrr=l WO M/s 
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FIGURE 17 

EQUI DENSITY PLOT (I0~14/m~3) 


CQ2; TEMP=300 K, V'ELOQTT=Q m/5 
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FIGURE 18 


EQUI DENSITY PLOT ( 1 1 5/m rt 3) 


C02; TEW P= 300 K, VELOCITY O M/5 
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FIGURE 19 


EQUIDENSITY PLOT ( 10~ 1 6/m~3) 


002; TEMP=300 K. VEL0CITt'=0 M/5 
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FIGURE 21 


EQUI DENSITY PLOT ( 1 CT 1 4/ 


C02; TEMP =10 K, VELOCITY^ 4000 M 
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Appendix A 


Nomenclature 


i 

$ 


6 

m 


k 

T 

M 

n=n(r,t) 


Qo(t) 


Cos 0 


t 

V r 

$ 

erfc 


= phase space particle density 
= position relative to source 
- velocity of particle in phase space 
= delta function 
= mass of molecule 
= Boltzmann's constant 
= Kelvin temperature 
= mean molecule velocity 
= number density 

= number flux rate from a source 

— i -4 
= ii*r 

= time source is on 
= del operator (coordinate) 

= velocity distribution function 
= complimentary error function 
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INTRODUCTION 

The self-induced molecular contamination around the Space Station could have 

, ... * ».<f ' 

adverse effects on Space Station components (for example solar panels) as well 
as scientific experiments (for example experiments involving photons) that 
might be done on or near the Space Station. One such source of self-induced 
contaminates is the propulsion system. Consideration is presently being given 
to continuously or semi -continuously firing small resistojets (in the hundreds 
of millinewtons thrust range) to overcome on-orbit frictional drag. 

This report outlines an initial effort to develop a mathematical model to 
study the effect of nozzle design on the induced molecular environment around 
the Space Station produced by simple gas propellants. The mathematical model 
would allow one to follow the expansion of the gas from the throat of a nozzle 
to the nozzle exit plane and then into the space external to the nozzle. A 
description of the proposed two-dimensional model and the mathematical 
calculations (by numerical methods) that could be performed in a rather 
complex computer code that must be written to obtain actual predictions are 
outlined in the remainder of this report. 

BASIC APPROACH 


Because of the expected relatively low densities of the propellant gas in the 
region of the exit plane of the nozzle, the fluid flow will be considered from 
the particle (molecular) nature of the gas using the methods.-iif^'/^ 
kinetic-statistical theory. The starting point for describing the gas flow 
will be the Boltzmann^ equation for binary collision. In particular, the 

Jt- 


model described below will be based on the steady state formulation of the BGK 

( 2 ) (31 

approximation' ' or "first iterate equation"' ' of the Boltzmann 
equation. Then by an iterative process a solution to the Boltzmann equation 
is obtained. 


MATHEMATICAL FORMULATIONS 


For the steady state the first iterate equation in differential form is 

V, f CC'i r) - /?^r) (f t ~f) 


(1) 
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where Maxwellian molecules have been assumed for describing the 
interparticle collisions. (See Appendix for equation symbol meanings.) With 
the appropriate boundary conditions, the first order solution to equation (1 ) 
can be obtained. Then the following first moments can be obtained: 


& f(r,£)J S 


* - m 


or~ 




3 

/ir- 


12 ) 


( 3 ) 


( 4 ) 


The value of^ , ^ , T, and -p can be the starting values for a second 
iteration, etc. 


Using the method of characteristics equation (1) for can be written in 
integral form as follows: - , . 

H?,*) -i 


( 5 ) 


The 


term is to be determined from one dimensional isentropic flow. 


SPECIFIC NOZZLE GEOMETRY AND BOUNDARY CONDITION ASSUMPTIONS 

The following section describes the suggested nozzle geometry and boundary 

conditions for^ which can be used to try to simulate the flow of a propellant 

through a nozzle as well as the flow of the plume. Figure 1 is a diagram of a 

two-dimensional nozzle and accompanying physical boundaries where -f-f must 

B 

be specified. 

The following boundary conditions {see Figure 1) for -f- seem appropriate: 


= £> for the inward component of velocity 


= £itT \\ e 
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fs, = ° 
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diffuse scattering 


where /r \ g and are nunfcer densities on the respective boundaries. In 

Figure 1 , the outer boundary Rj is assumed to be approaching infinity thus 
the number of particles with an inward component of velocity is zero. The 
number of particles on the outer surfaces (& 3 of this nozzle is assumed 

to be zero. The value of -f at is assumed to be that obtained from 
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isentropic flow conditions at the throat of the nozzle. In Figure 1 the 
inside nozzle wall boundary conditions are for complete diffuse scattering, 
but spectral conditions could also be applied. 

DISCUSSION 


Within the limitation of this model, one would be able to predict from throat 
conditions the number density at all points in space. Thus, with the model 
one can generate number densities at all points around the nozzle as well as 
number densities near the lip of the nozzle where back flow may be 





APPENDIX 
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AX 

u* 
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constant 

length in characteristic direction 

phase space particle density 

phase space particle density at boundary 

Maxwellian velocity distribution 

Boltzmann's constant 

mass of molecule 

number density 

nunber density at boundary 

To 

space position relative^nozzle 
location of boundary 
local temperature 
temperature at boundary 
temperature of gas at throat 
local mean velocity 
mean velocity at nozzle throat 
velocity 

number density at boundary 
del operator (coordinate) 
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